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ABSTRACT 

The discrepancy between the experimental decay rate of orthopositronium (o-Ps) 
and the QED theoretical prediction can be solved by invoking decays of o-Ps into 
exotic particles with branching ratios ~ 10 -3 . We show that considerations based 
on primordial nucleosynthesis and effective Lagrangians place a very stringent upper 
bound: B = T(o-Ps — > "exotic" + ...)/T(o-Ps) < 2 x 1(T 15 , ruling out the exotic 
decay solution to the puzzle. 



1. Introduction 

In the last few years, experiments |l], [| performed on the orthopositronium (o-Ps) system 
-e + e~ bound in the 3 Si state- have shown that its decay rate in vacuum is significantly 
higher than the QED prediction. The data obtained by Westbrook et al. |l|], give \ exp = 
7.0516 ± 0.0013 /is' 1 , while Nico et al. @ obtain X exp = 7.0482 ± 0.0016 us' 1 . However, the 
QED expectation || is X t h = 7.03830 ±0.00007 us' 1 , so that there is a significant discrepancy 
between theory and experiment. 

One of the proposed solutions to this problem is that, apart from the standard 37 decay, 
there are new disintegration channels of o-Ps. Indeed, if o-Ps has exotic decays, with a 
branching ratio on the order of 10 -3 , the theoretical prediction would increase in the right 
amount to be consistent with experiment. 

Much attention has been devoted to decays with final states that involve exotic particles. 
The decay mode o-Ps — > 7X, where X is a weakly interacting particle Q], or the decay into 
"invisible" final states, could solve the orthopositronium decay rate puzzle. Several groups 
have been searching for such new decays of o-Ps 0-0. The most restrictive experimental 
results on these exotic decays that have been published are: 

• o-Ps — > 7X has, if X is short-lived, the following upper limits on the branching ratios: 
2.0 x 10~ 4 for 847 keV < m x < 1013 keV §, 3.0 x 10" 4 for m x < 500 keV § 

• o-Ps — > 7X has, if X is long-lived, a branching ratio that is less than 1.1 x 10~ 6 , provided 
m x < 800 keV § 

• o-Ps — > "nothing" , has also been searched for. The absence of any "invisible" event 
gives a branching ratio less than 2.8 x 10 -6 ||. 

We see that the decay o-Ps — > "nothing" is ruled out as an explanation of the o-Ps 
anomaly, and that the decay o-Ps — > 7X is excluded in most of the possible X mass range, 
< nix < 1022 keV, but not in all the range. 

The experimental limit on o-Ps — > 7X is based on the search for a monochromatic photon 
in the final state. We note that a decay such as o-Ps — > 7X1X2, where there are now two 
exotics in the final state could provide a solution for the o-Ps puzzle, but there are no experi- 
mental limits on such decays, the reason being that now the photon is not monochromatic. 

In this letter we place bounds on the o-Ps exotic decays based on primordial nucleo- 
synthesis arguments. We will obtain much more restrictive bounds than the experimental 
results, for any decay of o-Ps into final states containing at least one exotic particle. 
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Simply stated, our idea is the following. A exotic particle that is produced in o-Ps decay 
is coupled to electrons, positrons and/or to photons with a strength that will determine 
the branching ratio for this mode. In the early Universe, that particle will be maintained 
in thermal equilibrium due to its exotic coupling to e + , e~ and/or 7. However, the exotic 
particle must decouple by the nucleosynthesis era (T ~ 1 MeV) since otherwise the successful 
predictions of the standard primordial nucleosynthesis would be spoiled. This requirement 
places strong limits on the exotic coupling, and consequently on the exotic decay of o-Ps. 

We are interested in the contribution to the effective degrees of freedom of new particles 
with masses up to the o-Ps mass, i.e., < m < 1022 keV. Therefore, in general, we cannot 
neglect the mass m compared to the temperatures T that we will have to consider, since the 
latter go down to T ~ 1 MeV. We will discuss how to treat this situation in Sect. 2. In Sect. 3 
we will discuss the 7X decay channel, with X a scalar particle, and in Sect. 4 we will discuss 
other decays. In Sect. 5 we will give our conclusions. 



2. The contribution of X to the effective degrees of freedom 

As we will see, we will be mainly interested in the exotic particle X being a scalar S. 
Let us then consider this case: X = S. The scalar particle S was in thermal equilibrium in 
the early Universe, due to its interaction with electrons, positrons and/or photons. It has an 
interaction rate T that can be calculated once this interaction is specified. 

As the Universe expands and cools, the interaction rate decreases and S decouples at a 
time to and at a temperature Tjj when 

T = H (1) 

where H is the expansion rate. Neglecting the curvature term, a completely justified approxi- 
mation in the early Universe, H is given by 

H > = ^ p = aAT) ^£ T < (2) 
3 H y K ' 90 K J 

In (fj) we have defined the effective degrees of freedom g*(T) (we follow the notation of Kolb 
and Turner 



As is well known, primordial nucleosynthesis offers a limit on the degrees of freedom 
contributing to the early Universe expansion at T ~ 1 MeV. It comes from a detailed 
comparison between the observed abundances of the primordial elements and the predictions 
of the Big Bang model in the nucleosynthesis era |TT[ . The limit is very stringent and places 



a bound on the energy density ps of S compared to the photon contribution p 7 [| 



Ps 



1/2 pi 



< 0.5 T ~ 1 MeV 



(3) 



The scalar particle 5 would of course contribute to g*. Its energy density is 

1 



Ps 



dpAnp 2 Ef(p 1 t) 



(4) 



where E = \Jp 2 + m 2 . Before decoupling, the particle is in thermal equilibrium and its 
distribution function / is that of a boson at temperature T 



f(p,t) = [exp(E/T)-l] 



-i 



T > T, 



D 



(5) 



The temperature T is a decreasing function of t. After decoupling, the number of S particles 
is conserved, which is expressed by 



f(p,t)=f(p D ,t D ) T<T D 

where the momentum has been redshifted as pR = pdRd, with R the scale factor. 
From (|5|) and (|6|) we get the expression 



(6) 



/CM) 



jp2 m 2 



D 



T <T, 



D 



(7) 



where we have defined the time dependent parameter 

= ^T D (8) 

The parameter 9 cannot be interpreted as the temperature of S, unless S is massless. 9 can 
be calculated using entropy conservation, which involves the function g*s{T), the effective 
degrees of freedom contributing to the entropy, 



g*s(T)R 3 T 3 = g* s (T D )RiT* 



3 t->3 



which gives 



(Here S refers to entropy). 



9 = T 



' g*s(T) 
^g*s{T D ) 



1/3 



(9) 



(10) 



1 A very recent analysis |12| on primordial nucleosynthesis tends to suggest that the limit is even smaller 
than ([|). Our final upper bound on exotic branching ratios would be then more stringent. 
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We have evaluated Ag*, as denned in (|J), numerically, substituting (0) and fliTf ) into (J|) 
and using the values of g*(T) and <?*s(T) calculated using their exact definitions (see for 



example [[HJ). Then, in order not to spoil the agreement of the predictions of the Hot Big 
Bang model with the observed primordial abundances, we impose the experimental limit (|3|). 
In the mass range < m < 1022 keV, the limit translates into the condition that the S 
species must decouple for all the temperatures in the range 

1 MeV < T < 100 MeV (11) 
3. The o-Ps — ► 7X decay, with X = S a scalar particle 

We will assume in this section that C and P are conserved. There are three ways in which 
this decay can then proceed: 

(A) S is coupled to electrons. (An equivalent model with X a pseudoscalar particle was 
considered by Samuel [|]].) The o-Ps decays via Fig. (l.A). The exotic coupling is described 
by a Yukawa-type interaction 

C a = 9aWS (12) 
where both here and in the following ip is the electron field. 

(B) S is coupled to photons. The decay proceeds through the diagram Fig. (l.B). The 
effective Lagrangian describing the new coupling is 

C B = \g B F^F^S (13) 

(C) There is a contact term coupling electrons to a photon and S as in Fig. (l.C). The 
Lagrangian describing the interactions has the gauge-invariant form 



Cc = \gc^^F^S (14) 



Let us start with model (A). We have an exotic particle S of mass 171$, with the interaction 
Lagrangian given in (|i~2"D. 

In the early Universe, the scalar S is in thermal equilibrium due to the processes 

7 5 -> e + e~ (15) 
eS — > e7 (16) 

where e can be either e + or e~. We have to consider all three processes since their contribution 
to the interaction rate are of approximately the same magnitude. 
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Let us work out the contribution of the process 7S — > e + e~ in some detail. We are 
interested in temperatures 3T > 3 To >> m e , ms- In this limit, the process has a total cross 
section given by 

2 

o-[ s ] = ^l ns / m 2 (17) 

The corresponding interaction rate (per unit time and per S) is 

r = n 7 < av > (18) 
where the number density of photons is given by 

n 7 = 2®T 3 (19) 



7T 



The thermalized quantity < av > appears in the evolution equation of the S number density 
ns (the superscript eq refers to the distribution in thermal equilibrium) 

h s + 3Hn s = - < av > n e J[n s - n e s g ] (20) 



One can find approximate expressions for < av > in |p^| . We write the expression for this 
quantity in the relativistic limit we are interested in 

< av >= f dn e y q dnf (1 - cos#) afiE^Esil - cos 9)} (21) 

Here Ei refers to the energy of the particle i and 9 is the angle between S and 7. Substituting 
now the cross section (|T7|) into this expression, one obtains after a straightforward calculation 

<->- ss^KM 1 -^ ^) +21n 3 (22) 

where £(x) is the Riemann £ function, £(se) its derivative and 7 is the Euler constant. 

We can now obtain restrictions on by requiring that S is decoupled in the temperature 
range determined in (|TT|). For each To in this range, we obtain an upper limit on g^, and 
our final result is necessarily the most stringent of them all. It turns out that this numerical 
limit is obtained at ~ 1 MeV. 

The formulae above actually only take account of the process (|I5D, which is one of the 
reactions that maintains S in equilibrium. The contributions of the other processes (|l^) can 
be worked out along the same lines. If one takes all these processes into account, the upper 
bound on the coupling gA appearing in the effective Lagrangian (O) is found to be 
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3A<5x 10- 10 = gT er (23) 
where we have defined (for a later purpose) the numerical value of the upper limit as g u ^ ver . 



With the Lagrangian (|T^), one can calculate the branching ratio 



_ r(oP-s -> jS) _ , 3 / ml \ 

Bs ~ T(oP-s) - 9a 8(tt* - 9)a* ^ ~ 4<) ( } 



This leads directely to one of our main results. Inserting (|23"D in (23), the following 
stringent limit is found 

B s < 2 x 1(T 15 f A (m s ) < 2 x 10~ 15 = B^ er (25) 

where the function fa (777,5) = Bg(ms) / 'Bs(ms = 0) = 1 — m s /Am 2 e < 1, with B s calculated 
using (pD. 



As we discussed at the beginning of this section, the decay o-Ps — > 75 is due to new inter- 
actions that can be described by effective Lagrangians. We have presented the consequences 
of assuming that the exotic decay is due to the new interaction (A). Using ( f[2f) we obtained 
the limit (p5|) on Bg. We have performed similar analyses in the cases (B) and (C), using 
(13) and (14) respectively. We will only quote the value of the limits. 



For the interactions (B), we obtain 

g B <2x 10- 7 GeV- 1 = gj per (26) 

and 

B s <2x 10- 17 f B (m s ) < 2 x 1(T 17 = B* pper (27) 
with f B (m s ) = (1 - ml/Amlf < 1. The limit is now obtained at T% ~ 100 MeV. 
For the case (C), we obtain 

5c <lx 10- 7 GeV" 2 = g u c pper (28) 

and 

B s < 1 x 10- 22 f c (m s ) ee B c apper (29) 
with f c {m s ) = (1 - m|/4m 2 ) 3 < 1. The limit is also obtained at T% ~ 100 MeV. 
The bound on Bs we finally can claim is the least restrictive, i.e., Bs < 2 x 10~ 15 . 
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4. Other o-Ps decays to exotic particles 

In order to discuss other decays, it will be very useful to understand why we have found 
that the limit on B$ is less restrictive in case (A) than in (B), and why the latter is less 
restrictive than case (C). 

The differences arise from the different dimensions D of the operators appearing in the 
effective Lagrangians. The operator in (|l2|) has D = 4, in fll3"l) it has D = 5, while in (jl4|) it 
has D = 6. Consequently, the coupling g& is dimensionless, while qb and gc have dimensions 
M~ l and M~ 2 , respectively. It is important to notice that there are two energy scales in 
our problem. On the one hand, we obtain restrictions at temperatures To, where decoupling 
is necessary in order not to excessively modify the predicted abundances in the primordial 
nucleosynthesis. On the other hand, we use these restrictions into the exotic o-Ps decay, that 
occurs at an energy scale E ~ m e . This leads to the following relations among the upper 
limits on the exotic couplings and on the branching ratios, that are valid up to an order of 
magnitude: 

B A ~ 1d±d B b B b ~ D D B c (30) 

upper 2 upper ' upper 2 upper \ J 

We can understand from this discussion that the least restrictive limit on B$ is in the case 
that the coupling is dimensionless, as it is in (A). This discussion will be the key to show the 
main point of this section, namely that all decays involving exotic particles in the final state 
have at least the upper limit B^ pper , corresponding to case (A) above. 

The o-Ps — > 7A decay, with X a pseudoscalar particle P, is for our purposes completely 
similar to the decay where X is a scalar particle. There are also three ways in which the 
decay can proceed, described by the three Lagrangians 

A = gxh^P (3i) 

£2 = \g2^sF^F^P (32) 

o 

£3 = \ 93^75^ F^P (33) 

The conclusions regarding the upper limits one can obtain on Bp are exactly the same as the 
corresponding cases (A), (B) and (C) above. Thus B£ r is the least restrictive upper limit 
we can place on the branching ratio of the decay o-Ps — ► 7P. One can think of interactions 
that do not conserve C and/or P and leading to o-Ps — > jX, with X a scalar or pseudoscalar 
particle. The results one finds are again exactly the upper limits that we found in cases (A), 
(B) and (C). 
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Let us now consider the decay o-Ps — > 7X, with X a spin one or two particle. The 
effective Lagrangian describing the exotic interaction contains operators with D > 4. The 
least restrictive limit will be obtained when D = 4, since then the coupling is dimensionless. 
In fact, since X now has more degrees of freedom that could contribute at the nucleosynthesis 
era, the decoupling temperature that we obtain is higher than the one we found in the case of 
only one degree of freedom (X — S or P) . Therefore the limit we obtain is even several times 
more restrictive than that of case (A). Thus B^ pper has to be seen as a conservative upper 
limit when X has spin different from zero. 

Another decay that has been experimentally investigated is o-Ps — > "nothing" . A possi- 
bility here is o-Ps — > SV, with both the scalar S and the vector V independently attached 
to electrons. Here we have two effective Lagrangians with D — 4, and two independent 
dimensionless couplings: 

£ 4 = g*WS (34) 
£5 = g^i^v, (35) 

Our arguments apply to the product of these two couplings. We have to take into account 
that one has now two exotic particles and thus more degrees of freedom, and the decoupling 
will have to occur at higher temperatures making the bound on the branching ratio T(o-Ps — > 
"nothing" )/r(o-Ps) even more stringent than B^ pper . Another possibility is to have contact 
terms like 

£ 6 = g^r^S (36) 
C 7 = g 7 ^a^(dMS (37) 

The coupling g e has dimension M" 1 and the coupling g 7 has dimension M~ 2 ; thus the bound 
on the branching ratio will be much more restrictive than B^ L . There are still other ways 
to obtain a o-Ps — > "nothing" decay. We will not list all the possibilities, since it is easy to 
see that we will always have limits more stringent than B^ pper . 

An interesting decay of o-Ps would be into a photon and two exotic particles. In fact, 
experiments are optimized to search for a monochromatic photon, so that they are not sensitive 
to this type of decay. We can apply our arguments based on nucleosynthesis to restrict 
this possibility. For example, consider o-Ps — > 7S , iS , 2 . The case that the two scalars are 
independently attached to electrons is similar to the above mentioned o-Ps — > SV decay, 
and one reaches the same conclusion, namely that B^ r is a conservative upper limit on the 
branching ratio of the exotic decay. If both scalar particles emerge in a contact term with 
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electrons, or if the scalars are coupled to two photons, the coupling is dimensionful and the 
bound is still more stringent. 

The type of arguments discussed in this section can be applied to any decay with exotic 
particles in the final state. One first has to write the effective Lagrangians that can lead to 
the decay channel. The order of magnitude of the limit on the branching ratio depends on the 
dimensions of the coupling, and to a lesser extent on the number of exotic degrees of freedom. 
The least stringent limit is obtained when the coupling is dimensionless and there is a scalar 
or pseudoscalar particle in the final state (model (A) or closely related), so that we place the 
following limit on any exotic decay 

„ r(o-Ps -> "exotic" + ...) i rr 

B = -i — — — — — '- < 2 x 1(T 15 38 

r o-Ps ~ v ; 



5. Conclusions 

The o-Ps decay puzzle would be solved if there were new decays with branching ratios 
~ 10~ 3 . Some of these decays have been ruled out by existing laboratory experiments. 

In this letter we have shown that all o-Ps decays involving exotic particles are restricted 
by nucleosynthesis arguments. Using an effective Lagrangian approach, we have been able to 
place a very stringent upper limit for any exotic decay 

B = r( °- Ps ^; WiC " + - ) < 2 x 10- 15 (39) 

r(o-Ps) ~ v ; 



Our result excludes the exotic decay solution to the o-Ps puzzle. 
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Figure Caption 

Fig.l. Exotic o-Ps decays due to the effective Lagrangians in models (A), (B) or (C). 



ff 



